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PREFACE 



Hyperbolic Partial Differential Equations (PDEs), and in particular first-order 
systems of conservation laws, have been a fashionable topic for over half a 
century. Many books have been written, but few of them deal with genuinely 
multidimensional hyperbolic problems: in this respect the most classical, though 
not so well-known, references are the books by Reiko Sakamoto, by Jacques 
Chazarain and Alain Piriou, and by Andrew Majda. Quoting Majda from his 
1984 book, “ the rigorous theory of multi-D conservation laws is a field in its 
infancy . We dare say it is still the case today. However, some advances have been 
made by various authors. To speak only of the stability of shock waves, we may 
think in particular of: Metivier and coworkers, who continued Majcla’s work in 
several interesting directions - weak shocks, lessening the regularity of the data, 
elucidation of the ‘block structure’ assumption in the case of characteristics with 
constant multiplicities (we shall speak here of constantly hyperbolic operators); 
Freistiihler, who extended Majcla’s approach to undercompressive shocks, of 
which an important example is given by phase boundaries in van der Waals fluids, 
as treated by Benzoni-Gavage; Coulombel and Secclri, who dealt very recently 
with neutrally stable discontinuities (2D- vortex sheets), thanks to Naslr-Moser 
techniques. 

Even though it does not pretend to cover the most recent results, this book 
aims at presenting a comprehensive view of the state-of-the-art, with particular 
emphasis on problems in which modern tools of analysis have proved useful. 
A large part of the book is indeed devoted to initial boundary value problems 
(IB VP), which can only be dealt with by using symbolic symmetrizers, and thus 
necessitate pseudo-differential calculus (for smooth, non-constant coefficients) 
or even para-differential calculus (for rough coefficients and therefore also non- 
linear problems) . In addition, the construction of symbolic symmetrizers conceals 
intriguing questions related to algebraic geometry, which were somewhat hidden 
in Kreiss’ original paper and in the book by Chazarain and Piriou. In this respect 
we propose here new insight, in connection with constant coefficient IBVPs. 
Furthermore, the analysis of (linear) IBVPs, which are important in themselves, 
enables us to prepare the way for the (non-linear) stability analysis of shock 
waves. In the matter of complexity, stability of shocks is the culminating topic 
in this book, which we hope will contribute to make more accessible some of 
the finest results currently known on multi-D conservation laws. Finally, quoting 
Constantin Dafermos from his 2000 book, 1 hyperbolic conservation laws and gas 
dynamics have been traveling hand-by-hand over the past one hundred and fifty 
years'. Therefore it is not a surprise that we devote a significant part of this book 
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to that specific and still important application. The idea of dealing with ‘real’ 
gases was inspired by the PhD thesis of Steplrane Jaouen after Sylvie Benzoni- 
Gavage was asked by his advisor, Pierre-Arnaud Raviart to act as a referee in 
the defense. 

This volume contains enough material for several graduate courses - which 
were actually taught by either one of the authors in the past few years ■ depend- 
ing on the topic one is willing to emphasize: hyperbolic Cauchy problem and 
IBVP, non-linear waves, or gas dynamics. It provides an extensive bibliography, 
including classical papers and very recent ones, both in PDE analysis and in 
applications (mainly to gas dynamics). From place to place, we have adopted an 
original approach compared to the existing literature, proposed new results, and 
filled gaps in proofs of important theorems. For some highly technical results, 
we have preferred to point out the main tools and ideas, together with precise 
references to original papers, rather than giving extended proofs. 

We hope that this book will fulfill the expectations of researchers in both 
hyperbolic PDEs and compressible fluid dynamics, while being accessible to 
beginners in those fields. We have tried our best to make it self-contained and 
to proceed as gradually as possible (at the price of some repetition), so that the 
reader should not be discouraged by her/his first reading. 

We warmly thank Jean-Frangois Coulombel, whose PhD thesis (under the 
supervision of Benzoni-Gavage and with the kind help of Guy Metivier) provided 
the energy necessary to complete the writing of the most technical parts, for his 
careful reading of the manuscript and numerous useful suggestions. We also thank 
our respective families for their patience and support. 



Lyon, April 2006 



Sylvie Benzoni-Gavage 
Denis Serre 
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INTRODUCTION 



Within the field of Partial Differential Equations (PDEs), the hyperbolic class is 
one of the most diversely applicable, mathematically interesting and technically 
difficult: these (certainly biased) qualifying terms may serve as milestones along 
an overview of the field, which we propose prior to entering the bulk of this book. 



Applicability. Hyperbolic PDEs arise as basic models in many applications, 
and especially in various branches of physics in which finite-speed propagation 
and/or conservation laws are involved. To quote a few, and nonetheless funda- 
mental examples, let us start with linear hyperbolic PDEs. The most ancient 
one is undoubtedly the wave equation - also known in one space dimension as 
the equation of vibrating strings - dating back to the work of d’Alembert in the 
eighteenth century, which is closely related to the transport equation. We also 
have in mind the Maxwell system of electromagnetism, as well as the equation 
associated with the Dirac operator. Theoretical physics is a source of several 
semilinear equations and systems - semilinearity being characterized by a linear 
principal part and non-linear terms in the subprincipal part - for example, the 
Klein/sine-Gordon equations, the Yang-Mills equations, the Maxwell system for 
polarized media, etc. The non-linear models - often quasilinear - are even more 
numerous. The most basic one is provided by the so-called Euler equations of gas 
dynamics, which opened the way (controversially) in the late nineteenth century 
to the shock waves theory (later revived, in the 1940s, by the atomic bomb 
research, and still of interest nowadays for more peaceful applications, in medicine 
for instance). Speaking of flows, a prototype of scalar, one-dimensional conser- 
vation law was introduced in the 1950s in traffic flow modelling (under some 
heuristic assumptions on the drivers’ behaviour), which is nowadays referred to 
as the Lighthill-Whitham-Richards model. Other non-linear hyperbolic models 
include: the equations of elastodynamics (of which a linear version is widely 
used, in the modelling of earthquakes as well as in engineering problems with 
small deformations); the equations of chemical separation (chromatography, 
electrophoresis); the magnetohydrodynamics (MHD) equations - the coupling 
between fluid dynamics and electromagnetism being quite relevant for planets 
and other astrophysical systems - the Einstein equations of general relativity; 
non-linear versions of the Maxwell system for strong fields, for example the 
Born-Infeld model. Hyperbolic equations may also arise as a byproduct of an 
elaborate piece of analysis, as in the modulation theory of integrable Hamiltonian 
PDEs (like the Korteweg-de Vries equation and some non-linear Schrodinger 
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equations), in which the envelopes of oscillating solutions are described by 
solutions of the (hopefully hyperbolic) Whitham equations. 

This list of hyperbolic PDEs is by no means exhaustive. Of course most of 
them are to some extent approximate: more realistic models should also involve 
dissipation processes (for instance in continuum mechanics) or higher-order 
phenomena, and thus be (at least partially) parabolic or dispersive. However, 
large-scale phenomena are usually governed by the hyperbolic part: the relevance 
of hyperbolic PDEs in many applications is in no doubt. 

Mathematical interest. For both mathematical reasons and physical rele- 
vance, hyperbolicity is associated with a space-time reference frame, in the sense 
that there exists a co-ordinate (most often the physical time) playing a special 
role compared to the other co-ordinates (usually spatial ones). Of course, changes 
of co-ordinates are always possible and we may speak of time-like co-ordinates 
and of space-like hypersurfaces: this terminology is familiar to people used to 
general relativity, and is also relevant in every situation where a hyperbolic 
operator is given. Except in one-dimensional frameworks, it is by no means 
possible to interchange the role of space and time variables: the distinction 
between time and space is a crucial feature of multidimensional hyperbolic PDEs, 
as we shall see in the analysis of Initial Boundary Value Problems. 

Multidimensional hyperbolic PDEs constrast with one-dimensional ones from 
several points of view, in particular in connection with the important notion of 
dispersion. Indeed, recall that the most visible feature of hyperbolic PDEs is 
finite-speed propagation. In several space dimensions, when the information is 
propagated not merely by pure transport, it gets dispersed: this dispersion of 
signals is itself responsible for a damping phenomenon in all L p norms with 
p > 2 (by contrast with what usually happens with the L 2 norm, independent of 
time by a conservation of energy principle), and is associated with special, space- 
time estimates called Strichartz estimates , obtained by fractional integration - 
Strichartz estimates have been proved much fruitful in particular in the analysis 
of semilinear hyperbolic Cauchy problems. 

Another point worth mentioning is the diversity of mathematical tools that 
have been found useful to the theory of (linear) multidimensional hyperbolic 
PDEs, ranging from microlocal analysis to algebraic topology (not to mention 
those that still need to be invented, as we shall suggest below!). The former has 
been widely used to study the propagation of singularities in wave-like equations. 
In the same spirit, pseudo- (or even para-) differential calculus is of great help 
to study linear hyperbolic problems with variable coefficients, as we shall see in 
the third and fourth parts of this book. The link to algebraic topology might 
seem less obvious to unaware readers and deserves a little explanation. When 
studying constant-coefficients hyperbolic operators we are led to consider, in the 
frequency space, algebraic manifolds called characteristic cones - which are by 
definition zero sets of symbols, and are linked to finite-speed propagation. The 
fundamental solution, say E, of a constant-coefficients hyperbolic PDE is indeed 
known to be supported by the convex hull of T, the forward part of the dual 
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of the characteristic cone. In some cases, it happens that E is supported by 
r only; the open set co(T)\r, on which E vanishes, is then called a lacuna. 
For example, the wave equation in dimension 1 + d with d odd and d > 3, has a 
lacuna: its fundamental solution is supported by the dual characteristic cone itself 
(this explains, for instance, the fact that light rays have no tail). The systematic 
study of lacunas is related to the topology of real algebraic sets. 

Compared to linear ones, non-linear problems display fascinating new fea- 
tures. In particular, several kinds of non-linear waves arise (shocks, rarefaction 
waves, as well as contact discontinuities). They are present already in one space 
dimension. The occurrence of shock waves is connected with a loss of regularity 
in the solutions in finite time, which can be roughly explained as follows: 
non-linearity implies that wave speeds depend on the state; therefore, a non- 
constant solution experiences a wave overtaking, which results in the creation of 
discontinuities in the derivatives of order m — 1, if m is the order of the system; 
such discontinuities are called shock waves, or simply shocks. After blow-up, that 
is after creation of shock(s), solutions cannot be smooth any longer. This yields 
many questions: what is the meaning of the PDEs for non-smootlr solutions; 
can we solve the system in terms of weak enough solutions, and if possible in a 
unique, physically relevant way? The answer to the first question has been given 
by the theory of distributions, which is somehow the mathematical counterpart of 
conservation principles in physics: conservation of mass, momentum and energy, 
for instance (or Ampere’s and Faraday’s laws in electromagnetism) make sense 
indeed as long as fields remain locally bounded. The drawback is - as has long 
been known - that weak solutions are by no means unique, and this seems to hurt 
the common belief that PDE models in physics describe deterministic processes. 
This apparent contradiction may be resolved by the use of a suitable entropy 
condition , most often reminiscent of the second principle of thermodynamics. 
In one space dimension, entropy conditions have been widely used in the last 
decades to prove global well-posedness results in the space of Bounded Vari- 
ations (BV) functions - a space known to be inappropriate in several space 
dimensions, because of the obstruction on the L p norms (see below for a few 
more details). Entropy conditions are expected to ensure also multidimensional 
well-posedness, even though we do not know yet what would be an appropriate 
space: one of the goals of this book is to present a starting point in this direc- 
tion, namely (local in time) well-posedness within classes of piecewise smooth 
solutions. 

Finally, the concept of time reversibility is quite intriguing in the framework 
of hyperbolic PDEs. On the one hand, as far as smooth solutions are concerned, 
many hyperbolic problems are time reversible, and this seems incompatible with 
the decay (already mentioned above) of L p norms for p > 2 in several space 
dimensions. This paradox was actually resolved by Brenner [22,23], who proved 
that multidimensional hyperbolic problems are ill-posed, in Hadamard’s sense, 
in L p for 2. Incidentally, Brenner’s result shows that the space BV , which is 
built upon the space of bounded measures, itself close to L 1 , cannot be appropri- 
ate for multidimensional problems. On the other hand, time reversibility is lost 
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(as a mathematical counterpart of the second principle of thermodynamics) once 
shocks develop, whence a loss of information, the backward problem becoming ill- 
posed. As a matter of fact, shocks may be viewed as free boundaries and as such 
they can be sought as solutions of (non-standard) hyperbolic Initial Boundary 
Value Problems (IB VP): it turns out that most of the well-posed hyperbolic 
IBVPs are irreversible, as will be made clear in particular in this book ■ a large 
part of this volume is indeed dedicated to a systematic study of IBVPs, either 
for themselves, or in view of applications to well-posedness in the presence of 
shock waves. 

Difficulty. Even when a functional framework is available, a rigorous analysis 
of hyperbolic problems often requires much more elaborate (or at least more 
technical) tools than for elliptic or parabolic problems, notably to cope with the 
lack of smoothing effects. The situation is even worse in the non-linear context, 
where functional analysis has been useless in the study of weak entropy solutions 
so far (except for first-order scalar equations). This is why our knowledge of 
global-in-time solutions is so poor, despite tremendous efforts by talented math- 
ematicians. Speaking only about the Cauchy problem for quasilinear systems of 
first-order conservation laws, in space dimension d with n scalar unknowns, we 
know about well-posedness only in the following cases. 

• Scalar problems (n = 1), thanks to Kruzklrov’s theory [105]. 

• One space dimension ( d = 1) and small data of bounded variation: existence 
results date back to Glimm’s seminal work [70]; uniqueness and continuous 
dependence have been obtained by Bressan and coworkers (see, for instance, 
[25-27]). 

• Small smooth data and large enough space dimension (for then dispersion 
can compete with non-linearity and prevent shock formation) : most results 
from this point of view have been established by Klainerman and coworkers. 
See, for instance, Hormander’s book [88]. 

Amazingly enough, none of these results apply to such basic systems as the full 
gas dynamics equations in one space dimension (n = 3, d = 1) or the isentropic 
gas dynamics equations (n = 2) in dimension d > 2. 

Other results solve only one part of the problem: 

• Global existence for general data when d = 1 and n = 2 (under a genuine 
non-linearity assumption) by means of compensated compactness. This was 
achieved by DiPerna [49], following an idea by Tartar [202], Solutions are 
then found in L°° . Unfortunately, no uniqueness proof in such a large space 
has been given so far, except for weak-strong uniqueness (uniqueness in L°° 
of a classical solution) . 

• Local existence of smooth solutions for smooth data. This is quite a 
good result since it shows at least local well-posedness. It is attributed 
to several people (Friedrichs, Garding, Kato, Leray, and possibly others), 
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depending on specific assumptions that were made. Unfortunately, its 
practical implications are limited by the smallness of the existence 
time - recall that shock formation precludes, in general, global existence 
results within smooth functions. 

Having this (modest) state-of-the-art in mind, we can foresee a compromise 
regarding multidimensional weak solutions and non-linear problems: it will con- 
sist of the analysis of piecewise smooth solutions (involving a finite number 
of singularities like shock waves, rarefaction waves or contact discontinuities), 
tractable by ‘classical’ tools. This is the point of view we have adopted here, which 
defines the scope of this book: we shall consider either (possibly weak) solutions 
of linear problems with smooth coefficients or piecewise smooth solutions of 
non-linear problems - Cauchy problems and also of Initial Boundary Value 
Problems - to multidimensional hyperbolic PDEs. We now present a more 
detailed description of the contents. 

We have chosen a presentation involving gradually increasing degrees of 
difficulty: this is the case for the ordering of the three main ‘theoretical’ parts 
of the book the first one being devoted to linear Cauchy problems, the second 
one to linear Initial Boundary Value Problems, and the third one to non-linear 
problems; this is also the case inside those parts - the first two parts starting 
with constant coefficients before going to variable coefficients, and the third one 
starting from Cauchy problems, then going to IBVPs, and culminating with the 
shock waves stability analysis. As a consequence, readers should be able to find 
the information they need without having to enter overcomplicated frameworks: 
most chapters are indeed (almost) self-contained (and as a drawback, the book 
is not free from repetitions). 

Another deliberate choice of ours has been to concentrate on first-order 
systems, even though we are very much aware that higher-order hyperbolic PDEs 
are also of great interest. This is mainly a matter of taste, because we come from 
the community of conservation laws. In addition, we think that the understanding 
of either one of those classes (first-order systems or lriglrer-order scalar equations) 
basically provides the understanding of the other class (see, for instance, the book 
by Chazarain and Piriou [31], Chapter VII). Consistently with that choice, the 
main application we have considered is the first-order system of Euler equations 
in gas dynamics, to which the fourth part of the book is entirely devoted. We 
have tried to temperate this ‘monomaniac’ attitude by referring from place to 
place to higher-order equations, and in particular to the wave equation, which is 
the source of several examples throughout the theoretical chapters. 

Finally, to keep the length of this book reasonable, we have decided not to 
speak of (nevertheless important) questions that are too far away from the shock 
waves theory. Thus the reader will not find anything about the propagation 
of singularities as developed by Egorov, Hormander and Taylor. Likewise, non- 
local boundary operators as they appear, for instance, in absorbing or trans- 
parent boundary conditions will not be considered, and all numerical aspects of 
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hyperbolic IBVPs will be omitted, despite their great theoretical and practical 
importance. 

First part. The theory of linear Cauchy problems is most classical, even 
though some results are not that well-known. The chapter on constant-coefficient 
problems is the occasion of pointing out important definitions: Friedrichs sym- 
metriz ability; directions of hyperbolicity; strict hyperbolicity and more generally 
what we call constant hyperbolicity - the eigenvalues of the symbol of a so- 
called constantly hyperbolic operator are semisimple and of constant multiplicity, 
instead of being simple in the case of strict hyperbolicity. Throughout the 
book, all hyperbolic operators will be assumed either Friedrichs symmetrizable 
or constantly hyperbolic (or both), as is the case for most operators coming 
from physics. The chapter on variable-coefficients Cauchy problems presents, 
in more generality, the symmetrizers technique, and in particular introduces 
the notion of symbolic symmetrizers, thus illustrating the power of pseudo- 
differential calculus (for infinitely smooth coefficients) and even para-differential 
calculus (for coefficients of limited regularity). 

Second part. The theory of Initial Boundary Value Problems (IB VP) is 
inspired from, but tremendously more complicated than, the theory of Cauchy 
problems. A kind of introductory chapter is devoted to the easier case of 
symmetric dissipative IBVPs. The second chapter addresses constant-coefficients 
IBVPs in a half-space, in which a central concept arises, namely the (uniform) 
Lopatinskii condition. This stability condition dates back to the 1970s: simultane- 
ously with a work by Lopatinskii ( [122], unnoticed in the West, Lopatinskii being 
more famous for his older work on elliptic boundary value problems [121]), it was 
worked out by Kreiss [103], and independently by Sakamoto [174] for higher-order 
equations; in acknowledgement of Kreiss’ work on first-order hyperbolic systems 
we shall rather call it the (uniform) Kreiss-Lopatinskii condition, and we shall 
also speak of Kreiss’ symmetrizers, which are symbolic symmetrizers adapted to 
IBVPs. The necessity of Kreiss’ symmetrizers shows up indeed when a Laplace- 
Fourier transform is applied to the equations (Laplace in the time direction 
and Fourier in the spatial boundary direction): to obtain an a priori estimate 
without loss of derivatives we need to multiply the equations by a suitable 
matrix-valued function, depending homogeneously on space-time frequencies - 
thus being a symbol - in place of the energy tensor of the symmetric dissipative 
case; that matrix-valued symbol is what we call a Kreiss symmetrizer. The 
actual construction of Kreiss’ symmetrizers is quite involved, and requires a 
good knowledge of linear algebra and real algebraic geometry. For this reason, 
a separate chapter is devoted to the construction of Kreiss’ symmetrizers. The 
interplay with algebraic geometry (formerly developed by Petrovskii, Oleinik 
and their school) is a deep reason why we need a structural assumption such 
as constant hyperbolicity: even with this, there remain tricky points to deal 
with, namely the so-called glancing points, where eigenvalues lack regularity. 
The chapter on variable-coefficient IBVPs focuses more on the calculus aspects 
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of the theory: it shows how to extend well-posedness results to more general 
situations - variable coefficients with either infinite or limited regularity, non- 
planar boundaries - by means of pseudo- or para-differential calculus. 

The remaining chapters of the second part are devoted to more peculiar 
topics: characteristic boundaries (which yield involved additional difficulties); 
homogeneous IBVPs (which turn out to require only a weakened version of the 
uniform Lopatinskii condition); the so-called class WR, which consists of certain 
< ^°°-well-posed problems and is generic in the sense that it is stable under small 
disturbances of the operators, but displays estimates with a loss of regularity. 
These topical chapters may be skipped by the reader insterested only in the 
applications to multidimensional shock stability. 

Third part. We must admit that the current knowledge of non-linear multi- 
dimensional hyperbolic problems is very much limited: all well-posedness results 
presented in this part are short-time results; nevertheless, their proofs are not 
that easy. A first chapter reviews Cauchy problems: symmetric (or Friedrichs- 
symmetrizable) ones, but also those with symbolic symmetrizers (at is the 
case for constantly hyperbolic systems), for which well-posedness was not much 
known up to now (the only reference we are aware of is a proceedings paper 
by Metivier [132]). Well-posedness is to be understood in Sobolev spaces of 
sufficiently high index, or to be more precise, in ff s (R d ) with s > d/2 + 1 (the 
condition ensuring that H s (M. d ) is an algebra, whose elements are at least 
continuously differentiable, by Sobolev’s theorem). In other words, we speak in 
that chapter only of smooth, or classical solutions, except in the very last section, 
where we recall the weak-strong uniqueness result of Dafermos and prepare the 
way for piecewise smooth solutions considered in the chapter on shock waves. 
Then ‘standard’ non-linear IBVPs are considered in a separate chapter, which 
is the occasion to see a simplified version of what is going on for shocks. The 
chapter on the persistence (or existence and stability) of single shock solutions 
was one of the main motivations to write this book. The idea was to give a 
comprehensive account of the work done by Majda in the 1980s [124-126], after 
it was revisited by Metivier and coworkers [56,131,133,134,136,140]. Initially, 
we intended to cover also non-classical (multidimensional) shocks, as considered 
by Freistiihler [58,59] and Coulombel [40]. But for clarity we have preferred to 
concentrate on Lax shocks, while avoiding as much as possible to use their specific 
properties so that interested readers could either guess what happens for non- 
classical shocks or refer more easily to [40] for instance. We have also deliberately 
omitted the most recent developments on characteristic and/or non-constantly 
hyperbolic problems. 

Fourth part. This concerns one of the most important applications of hyper- 
bolic PDEs: gas dynamics. In fact, the theory of hyperbolic conservation laws 
was developed, in particular by Peter Lax in the 1950s, by analogy with gas 
dynamics : terms like ‘entropy’, ‘compressive’ (or ‘undercompressive’) shock are 
reminiscent of this analogy, and the so-called Rankine-Hugoniot jump conditions 
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were initially derived (in the late nineteenth century) by these two engineers 
(Rankine and Hugoniot) in the framework of gas dynamics. There is a huge 
literature on gas dynamics, by engineers, by physicists and by mathematicians. 
In recent decades, the latter have had a marked preference for a familiar pressure 
law, usually referred to as the 7 -law, for it simplifies, to some extent (depending 
on the explicit value of 7 ), the analysis of the Euler equations of gas dynamics. 
We have chosen here to consider more general pressure laws, which apply to 
so-called real - at least more realistic - fluids and not only perfect gases (as was 
the case in earlier mathematical papers, by Weyl [218], Gilbarg [69], etc.). 

In a first chapter we address several basic questions, regarding lryperbolicity 
and symmetrizability. The second chapter is devoted to boundary conditions 
for real fluids, a very important topic for engineers, which has (surprisingly) not 
received much attention from mathematicians (see, however, the very nice review 
paper by Higdon [84]). 

This applied part culminates with the shock-waves analysis for real fluids, in 
the last chapter. Even though it seems to belong to ‘folklore’ in the shock-waves 
community, the complete investigation of the Kreiss-Lopatinskii condition for 
the Euler equations is hard to find in the literature: in particular, Majda gave 
the complete stability conditions in [126] but showed how to derive them only for 
isentropic gas dynamics; a complete, analytic proof was published only recently 
by Jenssen and Lyng [92]. By contrast, our approach is mostly algebraic, and 
works fine for full gas dynamics (of which the isentropic gas dynamics appear 
as a special, easier case). In addition, we give an explicit construction of Kreiss 
symmetrizers, which (to our knowledge) cannot be found elsewhere, and is fully 
elementary (compared to the sophisticated tools used for abstract systems). 

Fifth part. This is only a (huge) appendix, collecting useful tools and tech- 
niques. The main topics are the Laplace transform - including Paley-Wiener 
theorems - pseudo-differential calculus, and its refinement called para-differential 
calculus. Less space demanding (or more classical) tools are merely introduced 
in the Notations section below. 
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The set of matrices with n rows and p columns, with entries in a field K, is denoted 
by M, Jxp (K). If p — n, we simply write M„ (K). The latter is an algebra, whose 
neutral elements under addition and multiplication are denoted by 0 n and 
respectively. The space M nxp (IK) may be identified to the set of linear maps 
from K p to K". The transpose matrix is written M T . The group of invertible 
n x n matrices is GL„(K). If p = 1, M„ x i(K) is identified with K n . 

Given two matrices M, N G M„ (C), their commutator MN — NM is denoted 
by [M, N]. 

— -T — - 

If K = C, the adjoint matrix is written M*. It is equal to M , where M 
denotes the conjugate of M. We equip C m and R m with the canonical Hermitian 
norm 



ll*|| 




Or*a-) 1/2 . 



This norm is associated to the scalar product 



(x,y) = '^2, Xj y j = y*x. 
i 



The norm will sometimes be denoted |ar| , especially when a: is a space variable 
or a frequency vector (used in Fourier transform.) 

A complex square matrix U is unitary \iU*U = I n , or equivalently UU* = I n . 
The set U„ of unitary matrices is a compact subgroup of GL„ (C). Its intersection 
O n with M n (K) is the set of real orthogonal matrices. The special orthogonal 
group SO„ is the subgroup defined by the constraint det M = 1. 

As usual, M nxp (C) is equipped with the induced norm 



||M|| = sup 



\\Mx || 
11*11 



When the product makes sense, one knows that ||MAT|| < ||M|| ||JV||. When p = 
n, M„(C) is thus a normed algebra, and we have || M k \\ < ||Af|| fc . If Q is a unitary 
(for instance real orthogonal) matrix, one has ||Q|| = 1. More generally, the norm 
is unitary invariant, which means that ||M|| = ||PMQ|| whenever P and Q are 
unitary. 

If M G M„(C), the set of eigenvalues of M, denoted by SpM, is called the 
spectrum of M . The largest modulus of eigenvalues of M is called the spectral 
radius of M, and denoted by p(M). It is less than or equal to ||M||, and such 
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that 

p(M)= lim ||M fc || 1/fc . 

k — >-+oo 

The following formula holds, 

||M || 2 = p(M*M) = 

Several other norms on M„(C) are of great interest, among which is the Frobenius 
norm , defined by 



||M|| F := PT|m, fc | 2 . 

V i. fc 

Since ||M||| = Tr(M*M) = TV (MM*), we have ||M|| < ||M||f- 

A complex square matrix M is Hermitian if M* = M. It is skew-Hermitian if 
M* = —M. The Hermitian n x n matrices form an R-vector space that we denote 
by H„. The cone of positive-definite matrices in H„ is denoted by HPD„. When 
M is Hermitian, we have ||M|| = p(M). Every Hermitian matrix is diagonalizable 
with real eigenvalues, its normalized eigenvectors forming an orthonormal basis. 
The skew-Hermitian matrices with complex entries form an M- vector space that 
we denote by Skew„. We remark that M„(C) = H„ © Skew„ and Skew„ = 
iH n . The intersections of H„, HPD n and Skew„ with the subspace M„(R) of 
matrices with real entries are denoted by Sym n , SPD„ and Alt„, respectively. 
We have M„(R) = Sym ra © Alt„. Real symmetric matrices have real eigenvalues 
and are diagonalizable in an orthogonal basis. 

Given an n x n matrix M, one defines its exponential by 



exp M = e M 




which is a convergent series. The map 1 1 — > exp (tM) is the unique solution of the 
differential equation 



dA 
d t 



= MA, 



such that A(0) = It solves equivalently the ODE 



The exponential behaves well with respect to conjugation, that is 
exp (PMP- 1 ) = P(exp M)P ~ 1 

for all invertible matrix P. The eigenvalues of exp A are the exponentials of those 
of A. In particular, p(exp A) is the exponential of the maximal real part Re A, 
as A runs over SpA. The matrix exp(A + B) does not equal (exp A) (exp B) in 
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general, but it does when AB = BA. In particular, exp A is always invertible, 
with inverse exp(— A). Other useful formulae are 

exp(M T ) = (expM) T , exp M = exp M, exp(M*) = (exp M)*. 

The exponential of a Hermitian matrix is Hermitian, positive-definite. The map 

exp : H n — ► HPD„ 

is actually an homeomorphism. The exponential of a skew-Hermitian matrix is 
unitary. 

Let A G M n (C) be given. The space C n splits, in a unique way, as the 
direct sum of three invariant subspaces, namely the stable, unstable and central 
subspaces of A , denoted, respectively, E S (A), E U (A) and E C (A). Their invariance 
properties read 

AE S (A) = E s (A), AE u (A) = E u (A) and AE C (A) C E C (A). 

The stable invariant subspace is formed of vectors x such that (exp t A) x tends 
to zero as t — > +00, and then the decay is exponentially fast. The unstable 
subspace is formed of vectors x such that (exp tA)x tends to zero (exponentially 
fast) as t — > —00. The central subspace consists of vectors such that (expfH)a; is 
polynomially bounded on K. Since these spaces are invariant under A, this matrix 
operates on each one as an endomorphism, say A s , A U ,A C . The spectrum of A s 
(respectively, A u , A c ) has negative (respectively, positive, zero) real part. The 
union of these spectra is the whole spectrum of A, with the correct multiplicities. 
Hence the dimension of E S (A) is the number of eigenvalues of A of negative real 
part (these are called ‘stable eigenvalues’), counted with multiplicities. When 
E C (A) is trivial, meaning that there is no pure imaginary eigenvalue, A is called 
hyperbolic (in the sense of Dynamical Systems). 

Dunford— Taylor formula. Let 7 be a Jordan curve, oriented in the trigono- 
metric way, disjoint from SpH. Let a be the part of SpH that 7 enclose. Then 
the Cauchy integral 

defines a projector (that is P% = P a ) whose range and kernel are invariant under 
A. (Moreover, A commutes with P a ). The spectrum of the restriction of A to 
the range of P a is exactly the part of the spectrum of A that belongs to o. In 
other words, R(P a ) is the direct sum of the generalized eigenspaces associated 
to those eigenvalues in a. 

More information about matrices and norms may be found in [ 187 ]. 

Functional spaces 

Given an open subset f 2 of K 11 , the set of infinitely differentiable functions (with 
values in C) that are bounded as well as all their derivatives on O is denoted 
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by < ^’ b °°(r2). The set of compactly supported infinitely differentiable functions 
(also called test functions ) is denoted by @(f, l). Its dual !&'(£ 1) is the space of 
distributions. The derivation dj := d/dxj is a bounded linear operator on f2>(VL). 
Its adjoint is therefore bounded on The distributional derivative, still 

denoted by dj, is the adjoint of —dj. 

A multi-index a is a finite sequence (aq, . . . , a n ) of natural integers. Its length 
| a | is the sum JT a 3 . The operator 

d a := df 1 ■■■d“” 

is a derivation of order \a\. We also use the notation 

£Ot COt-n 

S — SI * S n 5 

when £ £ K". 

Given a function / : Q — > C, the differential of / at point X is the linear 
form 

n 

d f(X) : £ ~ d f{X)£ := £ djf(X). 

i = i 

The map A i— > d/(A) (that is the differential of /) is a differential form. The 
second differential, or Hessian of / at X is the bilinear form 



D 2 /(*) :(£,»?) ~ £ £* T]j didjf(X) . 

i,j = 1 

We may define differentials of higher orders D 3 /, . . . 

Given a Banach space IS, the Lebesgue space of measurable functions u : Q —> 
E whose pth power is integrable, is denoted by L P (Q; E). When E = R. or E = C, 
we simply denote L p (fl) if there is no ambiguity. The norm in L p (fl; E) is 



ml lp ■= 



t(a 



\\ P E dx ) 



^ 1/P 



If m £ N, the Sobolev space W m ’ p (fi; E) is the set of functions in L p (fl; E) whose 
distributional derivatives up to order m belong to L p . Its norm is defined by 

i/p 



\u\\w m <p := 



= E 



I LP 



0<*|<r 



If p = 2 and if IS is a Hilbert space, W m,2 (fl; E) is a Hilbert space and is denoted 
U m (fi; IS), or simply H m (Q) if IS = C or IS = R or if there is no ambiguity. 

Sobolev spaces of order s (instead of m) may be defined for every real 
number s. The simplest definition occurs when p = 2, fl = R" and E = C, where 
ff s (R") is isomorphic to a weighted space L 2 ((l + |£| 2 ) s d£) through the Fourier 
transform. For a crash course on H S (Q) (sometimes also denoted H S (Q)), we 
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refer the reader to Chapter II in [31]; for more details in more general situations, 
see for instance the classical monograph by Adams [1], The notation will 
stand for the Sobolev space H s equipped with the weak topology instead of the 
(strong) Hilbert topology. 

The Schwartz space of rapidly decreasing functions S fi (M n ) will simply be 
denoted by 5^ when no confusion can occur as concerns the space dimension. And 
similarly, its dual space, consisting of temperate distributions, will be denoted 
by S'". 

Other tools 

We have collected in the appendix various additional tools, ranging from standard 
calculus and Fourier-Laplace analysis to pseudo-differential and para-differential 
calculus: we hope it will be helpful to the reader. 
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LINEAR CAUCHY PROBLEM WITH 
CONSTANT COEFFICIENTS 



The general Cauchy problem 

Let d > 1 be the space dimension and x = (aq, . . . , Xd) denote the space variable, 
t being the time variable. The Cauchy problem that we consider in this section 
is posed in the whole space while t ranges on an interval, typically (0, T), 
where T < +oo. 

A constant-coefficient first-order system is determined by d + 1 matrices 
A 1 , . . . , A d , B given in M„(R), where n > 1 is the size of the system. Then the 
Cauchy problem consists in finding solutions u{x , t) of 



du 

dt 



d 

+ J2 Aa 

Oi—l 



du 

dx a 



Bu + /, 



(1.0.1) 



where f = f(x,t) and the initial datum u(-,t = 0)=a are given in suitable 
functional spaces. To shorten the notation, we shall rewrite equivalently 

d t u + ^ A a d a u = Bu + f. 

a 

When / = 0, the Cauchy problem is said to be homogeneous. A well-posedness 
property holds for the homogeneous problem when, given a in a functional 
space X, there exists one and only one solution u in Y), for some other 

functional space Y, the map 



X -> ^(0,T; Y) 



o^« 



being continuous. ‘Solution’ is understood here in the distributional sense. Exis- 
tence and continuity imply X C Y, since the map a i— > u(0) must be continuous. 
We use the general notation 

X^Y 

a i— > u(t). 

Since a homogeneous system is, at a formal level, an autonomous differen- 
tial equation with respect to time, we should like to have the semigroup 
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property 



S t+s = S t o S a , S,t> 0, 



this of course requires that Y = X. We then say that the homogeneous Cauchy 
problem defines a continuous semigroup if for every initial data a £ X, there 
exists a unique distributional solution of class < &’(IR + ; X). Note that the word 
‘continuous’ relies on the continuity with respect to time of the solution, but not 
on the continuity of 1 i— > St in the operator norm. Semigroup theory actually tells 
us that, if X is a Banach space, the continuity in the operator norm corresponds 
to ordinary differential equations, a context that does not apply in PDEs. 

When the homogeneous Cauchy problem defines a continuous semigroup on a 
functional space X, we expect to solve the non-homogeneous one using Duhamel’s 
formula : 

u(t) = S t a + [ S t - S f(s)ds, (1.0.2) 

Jo 

provided that at least / £ L 1 (0 ,T;X). For this reason, we focus on the homoge- 
neous Cauchy problem and content ourselves in constructing the semigroup. 

Before entering into the theory, let us remark that, since (1.0.1) writes 



du 

9t 



Pu + f , 



where P is a differential operator of order less than or equal to one, the 
order with respect to time of this evolution equation, the Cauchy-Kowalevska 
theory applies. For instance, if / = 0, analytic initial data yield unique analytic 
solutions. However, these solutions exist only on a short time interval (0,T*(a)). 
Since analytic data are unlikely in real life, and since local solutions are of little 
interest, we shall not concern ourselves with this result. 



1.1 Very weak well-posedness 

We first look at the necessary conditions for a very weak notion of well-posedness, 
where X = 5 ^(K d ) (the Schwartz class) and Y = the set of tempered 

distributions. Surprisingly, this analysis will provide us with a rather strong 
necessary condition, sometimes called weak hyperbolicity 1 . 

Let us assume that the homogeneous Cauchy problem is well-posed in this 
context. Let a be a datum and u be the solution. From the equation 



+ v A a — 
dt ^ dx a 
06—1 



= Bu, 



(1.1.3) 



^■Some authors call it simply hyperbolicity , and use the term strong hyperbolicity for the notion 
that we shall call hyperbolicity. Thus, depending on the authors, there is either the weak and normal 
hyperbolicities, or the normal and strong ones. 
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we obtain u £ < ^ oo (0 ,T\Y). This allows us to Fourier transform (1.1.3) in the 
spatial directions. We obtain that (1.1.3) is equivalent to 

r\ a d 

-^7 + i ^2 VaA a u = Bu. 

a=l 

Using the notation 

d 

Mv) ■= ^2v a A a , 

a — 1 

we rewrite this equation as an ODE in t , parametrized by 77 

^ = {B - iA(rj))u. 

Since u(-, 0) = a, the solution of (1.1.4) is explicitly given by 

u(r,,t) = e t(B - lA{ri »a( V ). 

By well-posedness (1.1.5) defines a tempered distribution for every choice of a in 
the Schwartz class, continuously in time. In other words, the bilinear map 

((/), ^>) 1 — > f , ip(r))*e t ( B ~ lA< '' n ^(l>(r)) dr), ( 1 . 1 . 6 ) 

Js. d 

which is well-defined for compactly supported smooth vector fields <f> and ip, is 
continuous in the Schwartz topology, uniformly for t in compact intervals. 

Let A be a simple eigenvalue of A(£) for some £ £ R d . Then, there is a < rf°° map 
(t, a) 1 — > ( /j,,r ), defined on a neighbourhood W of ( 0 ,£), such that n( 0 , £) = —i\ 
and 

( t 2 B - iA(a))r(t , a) = a)r(t, a ), ||r|| = 1. 

Let us choose a non-zero compactly supported smooth function 0 : — > C with 

0(0) 7 ^ 0. Then, for small enough t > 0, the condition 77 — £ Supp0 implies 

( t , t 2 ij) £ W. For such a t, we may define two compactly supported smooth vector 
fields by 

4? (v) ■= 0(V ~ t~ 2 Qr(t, t 2 if), ip\^) := 0(77 - t~ 2 £)£{t, f 2 ??), 

where £ is an eigenfield of the adjoint matrix (t 2 B — iA(a ))* , defined and 
normalized as above. We then apply (1.1.6) to (< 7 ^, t/? 4 ). The sequence (</> 4 ) t ^ 0 
is bounded in the Schwartz topology, and similarly is (^ 4 )i_> o- Therefore 

[ (^ 4 )*e 4 ( B - iA(T 'Vd77= f e' i( 4 ’ 42 ?))/ 4 (^-r)(t,t 2 77)|0(77-CA 2 )| 2 d77 

J R d JK. d 

is bounded as t — > 0. Since it behaves like cexp(— iX/t) for a non-zero constant c, 
we conclude that Im A < 0. Applying also this conclusion to the simple eigenvalue 
A, we find that A is real. 



(1.1.4) 

(1.1.5) 
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Linear Cauchy Problem with Constant Coefficients 



The case of an eigenvalue of constant multiplicity in some open set of 
frequencies p can be treated along the same ideas; it must be real too. Finally, the 
points p at which the multiplicities are not locally constant form an algebraic 
submanifold, thus a set of void interior. By continuity, the reality must hold 
everywhere. We have thus proved 

Proposition 1.1 The (d?’,S e ") well-posedness requires that the spectrum of 
A(£) be real for all f inR d . 

When (J^, 5C') well-posedness does not hold, a Hadamard instability occurs: 
for most (in the Baire sense) data a in and for all T > 0, the Cauchy problem 
does not admit any solution of class ^(0, T; =5^'). This is a consequence of the 
Principle of Uniform Boundedness. 

Example The Cauchy-Riemann equations provide the simplest system for 
which this instability holds. One has d = 1, n = 2: 

d t u i + d x u 2 = 0, d t u 2 - d x ui = 0. 



This example shows that a boundary value problem for a system of partial dif- 
ferential equations may be well-posed though the corresponding Cauchy problem 
is ill-posed. 

The converse of Proposition 1.1 does not hold in general, mainly because 
of the interaction between non-semisimple eigenvalues of A(£) with the mixing 
induced by B. Let us take again a simple example with d = 1, n = 2, and 



A = A 1 





Since the matrix 



exp(— if A) = I 2 — if A 



has polynomial growth, the Cauchy problem for the operator dt + Yla A a d a is 
well-posed in the sense, and even in the (J^, sense. Actually, its 

solution is explicitly given by 

ui(t) = ai — ta' 2 , u 2 (t) = a 2 . 



(We see that there is an immediate loss of regularity.) However, with our non-zero 
B, the matrix M := t(B — if A) satisfies M 2 = —it 2 fl 2 , which implies that 



exp (t(B — if A)) = cosw /2 + Sm W M, 

u> 



where u> = t{if ) 1 / 2 . Since 



Im ui = ±f 



1 1/2 




Strong well-posedness 
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we see that offdiagonal coefficients of exp M grow like exp(c|^| 1 / 2 ) as £ tends 
to infinity, provided t ^ 0. Then a calculation similar to the one in the proof 
of Proposition 1.1 shows that this Cauchy problem is ill-posed in the (dZ’jdZ') 
sense. 



1.2 Strong well-posedness 

The previous example suggests that the notion of well-posedness in the (rather 
weak) (S*, J7") sense might not be stable under small disturbance (the instability 
result would be the same with eB instead of B). For this reason, we shall merely 
consider the well-posedness when Y = X and X is a Banach space. We then 
speak about strong well-posedness in X (or A-well-posedness) . When this holds, 
the map St '■ a i— > u(t) defines a continuous semigroup on X. It can be shown 
that if X is a Banach space, there exist two constants c, u, such that 

\\S t \\c(x) < ce ut , (1.2.7) 

Proposition 1.2 Let X be a Banach space. Then well-posedness (with Y = X) 
for some B G M n (R) implies the same property for all B. 

This amounts to saying that well-posedness is a property of (A 1 , . . . , A d ) alone. 

Proof Assume strong well-posedness for a given matrix B 0 . The problem 



y- a a 9u 
dt dx a 

Ot—1 



BqU 



( 1 . 2 . 8 ) 



defines a continuous semigroup (S t )t>o- One has (1.2.7) with suitable constants 
c and u). From Duhamel’s formula, (1.1.3) with a matrix B = Bq + C instead of 
Bo, is equivalent to 

u(t) = St.a+ f St~ s Cu(s)ds. (1.2.9) 

Jo 

Then we can solve (1.2.9) by a Picard iteration. Let us denote by R.u the right- 
hand side of (1.2.9), and / = (0,T) (with T > 0) a time interval where we look 
for a solution. Because of (1.2.7), there exists a large enough N so that R N is 
contractant on ^(/;A). Therefore, there exists a unique solution of (1.1.3) in 
^(/; X). Since T is arbitrary, the solution is global in time. □ 



1.2.1 Hyperbolicity 

We first consider spaces X where the Fourier transform defines an isomorphism 
onto some other Banach space Z. Typically, X will be a Sobolev space H s (R d ) n 
and Z is a weighted L 2 -space: 

Z = L 2 s (W 1 ) n , L 2 (R d ) := {v € L( oc (R d ) ; (1 + |^| 2 ) s/2 u G L 2 (R d )}. 




